RF bifurcation of a Josephson junction : microwave embedding circuit requirements 



O 
O 

(N 

O 

Q 

(N 
(N 

O 

o 

di- 
ss 

c3 



I 

C 

o 
o 



> 
in 

(N 

o 

-I— > 



I 

C 

o 
o 



X 



V.E. Manucharyan, E. Boaknin, M. Metcalfe, R. Vijay, I. Siddiqi0 and M. Devoret 
Department of Applied Physics, Yale University, New Haven, CT 06511, USA 

(Dated: February 4, 2008) 

A Josephson tunnel junction which is RF-driven near a dynamical bifurcation point can amplify 
quantum signals. The bifurcation point will exist robustly only if the electrodynamic environment 
of the junction meets certain criteria. In this article we develop a general formalism for dealing with 
the non-linear dynamics of Josephson junction embedded in an arbitrary microwave circuit. We find 
sufficient conditions for the existence of the bifurcation regime: a) the embedding impedance of the 
junction need to present a resonance at a particular frequency lor, with the quality factor Q of the 
resonance and the participation ratio p of the junction satisfying Qp 2> 1, b) the drive frequency 
should be low frequency detuned away from lor by more than ^/Zlor/ {2Q). 

PACS numbers: 85.25.Cp, 84.30.Le, 84.40.Az, 84.40.Dc 



I. INTRODUCTION 

Amplifying very small electrical signals is an ubiqui- 
tous task in experimental physics. In particular, cryo- 
genic amplifiers working in the microwave domain found 
recently a growing number of applications in mesoscopic 
physics, astrophysics and particle detector physics^. We 
have recently proposed^ to use the dynamical bifurca- 
tion of a RF-biased Josephson junction as a basis for 
the amplification of quantum signals. A bifurcation phe- 
nomenon offers the advantage of displaying a diverging 
susceptibility which can be exploited to maximize the 
amplifier gain without necessarily sacrifying its band- 
width. Among all very low noise and fast solid-state 
microwave devices, the Josephson junction distinguishes 
itself by offering strongest non-linearity combined with 
weakest dissipation. However, these characteristics are 
not by themselves sufficient. The electrodynamic envi- 
ronment of the junction must also satisfy a certain num- 
ber of conditions in order for a controllable and minimally 
noisy operation to be possible. In the recent Joseph- 
son bifurcation amplifier experiments^, the junction was 
shunted by a lumped element capacitor. A large capac- 
itance had to be fabricated very close to the junction 
to minimize parasitic circuit elements, at the cost of se- 
vere complexity of patterning and thin-film deposition. It 
would be very beneficial experimentally to simply embed 
the Josephson junction in a planar superconducting mi- 
crowave resonator. The aim of this article is to establish 
theoretically the requirements that need to be imposed 
on the embedding impedance of the junction in order to 
obtain a bifurcation whose characteristics are suitable for 
amplification. 

The article is organized as follows: after having briefly 
indicated the connection between a bifurcating dynami- 
cal system and amplication, we review the simplest non- 
linear dynamical system exhibiting the type of bifurca- 
tion we exploit, namely the Duffing oscillator. We then 
describe the parameter space of the oscillator, focussing 
on the neighborhood of the first bifurcation and dis- 
cussing why this is the most useful region. Having laid 
the general framework for the analysis of our problem, 



we then consider the simplest practical electrical imple- 
mentation of the Duffing oscillator, a Josephson junction 
biased by an RF source through an arbitrary microwave 
circuit. The notion of embedding impedance is intro- 
duced. For concreteness, we first examine the particu- 
lar cases where the embedding impedance corresponds to 
simple series or parallel LCR circuits. This allows us to 
formulate the conditions under which the resulting non- 
linear electrical system can be mapped into the Duffing 
model. We then examine the arbitrary impedance case, 
finding that it must correspond to that of a resonator 
with an adequate quality factor. We end the article 
by discussing possible detailed experimental implemen- 
tations of resonators and a concluding summary. 



II. AMPLIFYING WITH THE BIFURCATION 
OF A DRIVEN DYNAMICAL SYSTEM 



Amplification using a laser, a maser or a transistor 
is based on energizing many microscopic systems, like 
atoms in a cavity or conduction electrons in a channel, 
each one being weakly coupled to the input signal. The 
overall power gain of the system, which is determined 
by the product of the number of active microscopic sys- 
tems and their individual response to the input parame- 
ter, can be quite substantial. However, noise can result 
from the lack of control of each individual microscopic 
system. This article explores another strategy for am- 
plification which involves a single system with only one 
very well controlled collective degree of freedom, which is 
driven to a high level of excitation. Here, the input signal 
is coupled parametrically to this system and influences its 
dynamics. The best known device exploiting this strat- 
egy is the SQUIDS but other devices of the same type 
have been proposed 7 . Let us discuss the general question 
of the gain (ratio of output to input) in such a system. 

A driven dynamical system such as a SQUID is gov- 
erned by a force equation which, quite generally, can be 
written as 
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T (Jt, Y, X, a) = F ext (t) (1) 

where X is the system coordinate, F ext (t) is a peri- 
odic external drive pumping energy in the system, a a 
parameter of the system and T a function describing its 
dynamics which is necessarily dissipative since informa- 
tion is flowing away to the next stage of amplification. 
We are interested in steady state solution of Eq. (P) , in 
which the energy flowing in from the source is balanced 
by the energy losses. In the example of the rf SQUID, 
X is the total flux through the SQUID loop, a the signal 
flux, and F the external driving flux with a frequency 
in the MHz range. For the dc SQUID 9 , the frequency 
of the external drive current is and X is the common 
mode phase difference while a is the flux through the loop 
formed by the junctions. In this article we also consider a 
Josephson-junction-based device like a SQUID, but it is 
driven by a rf signal at microwave frequencies to increase 
speed and does not have intrinsic dissipation. 

When we use the dynamical system as an amplifier, 
we are linking the input and output signals to the pa- 
rameters a and variable X, respectively. Specifically, the 
signal s(t) at the input of the amplifier induces a varia- 
tion Sa — As. For a small input, the output S (t) of the 
amplifier will depend linearly on the modification of X: 
S (t) = £{X(t) — X(t) s= o}. Since we are looking for a 
maximal signal gain S/ s, it is natural to find an oper- 
ating point where a small change in the parameter a is 
going to induce a large change in the dynamics of the sys- 
tem, provided we can keep all other parameters constant. 
The largest susceptibility is found at a saddle-node bifur- 
cation point and it is in the neighborhood of such points 
that we will operate the amplifier. The saddle-node bifur- 
cation occurs when the drive parameters exceed certain 
critical values. Previously, Yurke et al& have studied 
Josephson systems mostly in the regime beneath these 
critical values. Here, we consider similar systems, but 
we exploit instead the bistable regime beyond the critical 
values and the large susceptibilities accompanying it. In 
the next section, we examine a simple model exhibiting 
such a saddle-node bifurcation phenomenon. 



III. ONE MINIMAL MODEL FOR A 
BIFURCATING NON-LINEAR SYSTEM: THE 
DUFFING OSCILLATOR 

One of the most minimal model displaying the bifurca- 
tion phenomenon needed for amplification is a damped, 
driven mechanical oscillator with a restoring force dis- 
playing cubic non-linearity. The equation of this model, 
often called the Duffing linear+cubic oscillato r 10 ! 11 ' 12 , is 

mX + 7 Y + kX(l - vX 2 ) = Fcosuj d t + F N (t) (2) 



where X is the position coordinate of the mechanical 
degree of freedom, m its mass, 7 its damping constant, 
k the stiffness constant of the restoring force, v the non- 
linearity parameter. The right handside parameters F 
and u> are the amplitude and angular frequency of the 
driving force, respectively. For completeness, we have 
also added on the right handside a noise force term Fn (t) 
whose presence is imposed by the fluctuation-dissipation 
theorem. It defines, through its correlation function, a 
thermal energy scale for the problem, but in the follow- 
ing sections, we are going to assume that this scale can 
be made much smaller than all the other scales in the 
problem. The effect of fluctuations will be treated in a 
latter article. 

This simplification made, we can rescale the problem, 
changing the position and time coordinates, obtaining in 
the end a three-parameter model 

x 

x + — + x(l - x 2 ) — /cos wt (3) 

where the dot now refers to differentiation with respect 
to the rescaled time r. The following equations express 
the rescaled quantities in terms of the original ones: loq = 

-±=, T = LU t, w = Q = X = JDX, f = 

m 

For reasons which will become clear later, we want to 
consider the low damping limit of such a model and we 
suppose that the quality factor satisfies Q 3> 1. 

In the weak non-linear regime, i.e. x (t) <C 1, the 

frequency wo of natural oscillations (/ = 0) decreases 

• 2 2 4 

with the oscillation energy u = ^- + ^ — x as 

-o = l~y + 0( U 2 ) (4) 

This corresponds to the softening of the restoring force 
as the amplitude of the oscillatory motion increases. 

The bifurcation phenomenon can be described crudely 
as follows: if we drive the system below the small oscilla- 
tion resonance frequency, increasing the drive amplitude 
slowly, the resulting oscillation will be very small at first. 
However, at a certain drive strength, the system becomes 
unstable and tends to switch to a high amplitude oscil- 
lation where it can better meet the resonance condition. 

In a quantitative treatment of the weak non-linear 
regime, we seek a solution involving only the first har- 
monic of the drive frequency^ 

x (r) = ^x (r) e tWT + c.c. (5) 

where the time dependence of the complex harmonic 
amplitude x (r) is slow on the time scale w~ 1 . Retain- 
ing in the equation only the terms evolving in time like 
e™ T and ignoring d 2 x (r) /dr 2 , one finds the following re- 
lation for x(t)1£: 
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FIG. 1: Upper panels: Response of the Duffing linear+cubic oscillator (see Eq. [2]) as a function of the dimensionless detuning 
fi, for different values of the dimensionless driving amplitude /. The three panels correspond to three different values of the 
small oscillation quality factor Q = 20, 50, and 200 and is shown over the same reduced frequency range 0.9 < w < 1.05. The 
response for the critical drive amplitude f c , where the curve presents a single point with diverging susceptibilities dx/dQ. and 
dx/df, is shown in green. From bottom to top, the increments in the drive amplitude / are 2 dB (left), 2 dB (center) and 
3 dB (right). In the grey region, the oscillation amplitude is taking place in the strongly non-linear regime (see text) and the 
curves, which are calculated within the weak non-linear hypothesis, are not to be trusted. Lower panels: Stability diagram for 
the dynamical states corresponding to the system response shown in the upper panels. The y-axis is the drive amplitude /, 
scaled by the critical amplitude f c . The blue (/s(f2)) and red lines delimit the region of bistability for the system and 

correspond to the points of diverging susceptibility which are visible in the upper panels. The two curves meet at the critical 
point shown by a green dot, whose position is determined by the critical detuning f2 c = — v3 and drive amplitude f c . The 
dashed line (/ m ,(!l)) "continues" the blue and red lines and corresponds to the points of maximal susceptibility with respect 
to driving force. The grey regions correspond to that in the upper panels. Note that the "trusworthy" domain of bistability 
increases monotonously with the quality factor Q. 
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in which we have introduced the reduced detuning Q 

= 2Q(w- 1) 

The static solutions (x = 0) for the modulus x = |x| of 
the fundamental amplitude can be obtained as a function 
of the parameters (CI, /) for a given Q by solving the 
equation 
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The susceptibility is given by the implicit expression: 
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In the upper panels of Fig. [TJ we show x as a function 
of O for increasing values of / and for Q = 20, 50, 200. 
For small drive, the curve is the familiar Lorentzian re- 
sponse of an harmonic oscillator, displaying a maximum 
response on resonance at f2 = and a half-width at 
half- maximum (HWHM) point at f2 = 1. As the drive 
strength is increased, the resonance curve bends towards 
lower frequencies, an indirect manifestation of Eq. 
There is a critical drive f c at which appears for the first 
time a critical reduced detuning Q c such that the suscep- 
tibility dx/df diverges^ 2 -. We call x c the oscillation am- 
plitude at this critical point. Analytic calculations^^ 
lead to 
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To be consistent with our weak non-linear regime hy- 
pothesis, we must have x c -C 1 which implies in turn 
Q > 1. 

For drives / > f c , the response curve x (CI) develops an 
overhanging part in which there are three possible values 
for x at each value of CI. The smallest and highest values 
correspond to two metastable states with different oscil- 
lation amplitudes, whereas the intermediate value corre- 
spond to an unstable state for the system. We denote 
as fg (Cl) and fg (Cl) the boundaries of this bistability 
interval: fg is the force at which the system, submitted 
to an increasing drive with a fixed frequency will switch 
from the low to the high amplitude state. Starting from 
this state and decreasing the amplitude of the oscilla- 
tory force, the system will switch back to the low ampli- 
tude state at fg. This possibility of the Duffing system 
to "bifurcate" between two different dynamical states at 
fg (Cl) and fg (Cl) is the phenomenon we are exploiting 
for amplification and whose electrical implementation is 
the main topic of this paper. It is easy to see that any 
input parameter coupled to k or m in Eq. @ will induce 
variations of the line fg (Cl). Fixing the drive parameters 
in the vicinity of this line, very small changes in the input 
will induce large variations in the oscillation amplitude. 
The variations can be reversible if we chose a point to 
the right of the critical point (continuous amplifier oper- 
ation) or the variations can be hysteretic if we chose a 
point to the left of the critical point (latched threshold 
detector operation). 

In the limit Q > 1, analytic calculations can be carried 
further and lead toi 3 . 
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We define f rns (Cl) as the line of maximum susceptibility 
<9x/<9/ on the low- frequency side of the resonance curve. 
It defines the line of highest amplification gain below the 
bifurcation regime^. Its expression is given by: 



f ma (Cl) 3V3 Cl'/ 2 
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The susceptibility on the high frequency side of the 
critical point diverges as: 



<9x 
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fg(Cl) and f rns in the parameter space (fi, f / f c ) are inde- 
pendent of the parameters of Eq. |(3J) and can be deemed 
" universal" . 

The dynamical critical point (Cl = Cl c ,f/f c = 1) is 
found at the junction between the dashed line and the two 
bifurcation lines. One can develop an analogy between 
the parameter space (CI, f / f c ) and the phase diagram 
of a fluid undergoing a liquid-vapor transition, the dy- 
namical critical point corresponding to the critical point 
beyond which vapor and liquid cannot be distinguished 
by a transition (supercritical fluid regime), and the bifur- 
cation lines corresponding to the limit of stability of the 
supercooled vapor and superheated fluid on cither side 
of the 1st order transition line (spinodal decomposition 
phenomenon). 



A. Weak and strong nonlinear regimes for the 
simple Duffing equation 

Let us now further discuss the small amplitude condi- 
tion x <C 1, which is necessary for the above results to 
hold. In Appendix A we show that as long as 



— ^x 2 < 1; 1/2 < w < 1 



(10) 



the Duffing model has stationary solution of the form 



x{T) = \Y,[x 2k . l e^ k -^ T + c.c. 

k=l 



with the coefficients X2k-i decreasing with the order k. 
Only odd multiple of the drive frequency thus appear 
in this series. The first harmonic coefficient x\ is given 
by the stable solution of Eq. ([6]) in the limit |ix| <^ 1- 
Inequalities () 1 0(1 defines rigourously the weak non-linear 
regime. 

By contrast, in the strong non-linear regime x 2 > 
1, even harmonics start to proliferate as the oscillation 
amplitude increases, leading eventually to chaotic behav- 
io r 15 i 16 ' 17 . It is important to note that the SQUID does 
not avoid this regime, even if, in general, strong dissipa- 
tion prevents fully developed chaos in this device. 

Wanting at all cost to minimize noise in our use of 
this dynamical system for amplification, we want to avoid 
the strong non-linear regime. Keeping in mind that we 
are going to work with a small detuning (ui — luq) /u)q, 
a conservative boundary separating the weak from the 
strong non-linear regime can be introduced in parameter 
space by requiring 



where ACI is defined by Cl = Cl c — Af2 and tt c A57. 

In the lower panels of Fig. [T] we plot the bifurcation 
forces fg (blue) and fg (red), and f ms (dashed) normal- 
ized to the critical force as a function of the reduced drive 
frequency O. Note that the lines representing /g(O), 



x < 0.5 



(11) 



In the lower panels of Fig. [TJ the grey region corre- 
sponds to condition (JTTJ) being violated for at least one 
of the oscillation states. A lighter shade of grey marks 
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FIG. 2: (a) Schematic of a Josephson junction RF-biased 
though an arbitrary electric quadrupole containing only lin- 
ear components. The Thevenin (b) and Norton (c) represen- 
tations of the circuit of (a) shows the dipole seen by the junc- 
tion. Note that the new source amplitude may now depend 
on frequency. 



the hysteretic region between and /s(^) to indi- 

cate that the low amplitude state does not violate (fTTj) 
while the high amplitude does. 

Note that in the reduced parameter space (f2,/// c ), 
the line / S (Q) corresponding to condition (fTTj) has a 
rather drastic dependence on Q, in contrast with the 
other lines. Of course, if we would plot the stability 
boundary lines in the absolute parameter space (u,F), 
the line corresponding to would be fixed while the 
critical point (u) c , F c ) would strongly depend on Q. 

Whatever the representation, the important message 
which arises from the stability diagram is that the 
amount of "real estate" in parameter space that can be 
used for bifurcation amplification increases with Q. All 
points along the line located between f c and f s (£i) 

are potentially useful. By realizing a large enough Q, one 
can always "buy" the necessary amount of "real estate" 
in the stability diagram, irrespectively of the value of the 
other parameters of the system. Of course, higher Q will 
tend to lower the bandwidth, but we can compensate this 
effect by increasing the operating frequency of the device. 



and applying either the Norton or Thevenin representa- 
tion, we arrive at the circuits of Fig. and c respec- 
tively, where Y e (ui) and Z e (u>) are the admittance and 
impedance of the quadrupole seen from the junction. In 
this transformation of the problem, the amplitude of the 
source may now depend on frequency, but we ignore this 
complication for simplicity. It turns out that in the cases 
of interest, either the Norton or Thevenin (but not both) 
representions will satisfy this hypothesis adequately. 

The Josephson electrical element, represented by a 
cross in our schematics, is defined by its constitutive 
equation 



I(t)=I„ sin (*j(*)M) (12) 

involving the generalized flux defined by $j = 
/_ V(t')di^ and where (po — h/2e is the reduced flux 
quantum. Here V(t) and I(t) are respectively the voltage 
across the junction and current through it. 

Comparing Eq. (| 1 2[) with the constitutive equation of 
an inductance 



di 



understand 
-l 



J(*) = ~<f(t) 



why the quantity Lj = 
— <po/Io is referred to as the effec- 



*j=0, 

tive Josephson inductance. 

Note that Eq. (fT2Tj provides a non-linear link between 
$j and I and is, in some abstract sense, analogous to 
the relation between force and position for a non-linear 
restoring force. However, the external circuit described 
by Y e (ui) or Z e (u>) also participate in the restoring force 
and we have to go through one further step in order to 
establish a link between our electrical system and the 
Duffing oscillator model. 



A. Separation of the linear and nonlinear 
contributions of the Josephson element 



IV. THE RF-BIASED JOSEPHSON JUNCTION 

We now apply these general considerations to the prac- 
tical case of a Josephson junction (JJ) biased by an RF 
source. At RF frequencies, the impedance of the bias- 
ing circuit cannot be taken either as zero (ideal voltage 
source) or infinite (ideal current source). In general, as 
depicted by Fig. [5^, the biasing circuitry should be mod- 
eled as a specific, frequency-dependent linear quadrupole 
connecting the Josephson element of the junction to an 
ideal voltage source (the other circuit element of the 
junction, its linear capacitance, has been lumped in this 
quadrupole). Taking the point of view of the junction 



It is useful to split the Josephson element into its 
purely linear and nonlinear components. The linear 
contribution of the Josephson element is the impedance 
iLjuj — Zj(uj), and can be incorporated in the biasing 
impedance. The non-linear contribution, however, can 
be defined only by first referring to either the Norton or 
the Thevenin representation. 

We thus expand Eq. (fT2"|) in two different ways 

I(t) = i-# (t) - $j (tf + O [&(t)} (13) 

and 
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FIG. 3: Two possible ways of separating the linear and non- 
linear contributions of a Josephson junction to a microwave 
circuit. Different symbols were chosen for the nonlinear el- 
ements to suggest their different current-voltage relations. 
The Josephson inductance in both cases has the usual value 
Lj = (fio/Io- 



<f>j(t) = LjI(t) + ^I S (t) 



O [I 5 (t)} (14) 



corresponding to these two representations, respec- 
tively. Relations (fl"3]) and (fT4|) correspond respectively 
to a parallel and series combinations of a usual linear in- 
ductance of value Lj and a nonlinear element which is 
defined by the higher order terms in the equations (see 
Fig [3]) and which is also characterized by the parameter 
Lj. We will call these new components parallel nonlin- 
ear element (PNL) and series nonlinear element (SNL) 
respectively. They are represented by spider-like symbols 
in Fig. [3 

For the purpose of this paper, it will be sufficient to 
keep only the first nonlinear term in each of the expan- 
sions above. The cut-line spider symbol corresponds to 
the constitutive equation 




X 



(b) 
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FIG. 4: Equivalent Norton (a) and Thevenin (b) represen- 
tations of a linear circuit driving a PNL and SNL Josephson 
element respectively. The appropriate bias source is shown as 
either a parallel current source or a series voltage source. 




lit) 



(15) 



while the spine- line spider symbol corresponds to the con- 
stitutive equation 



FIG. 5: Josephson junction biased by (a) a parallel LCR 
circuit and (b) a series LCR circuit. In the first case, 
L~ 1 = Lpl + Lj 1 and in the second case, L s — L se + Lj 
where L pe and L 3e are the inductances contributed by the en- 
vironment. 



*(*) = S 4/ 3 W (16) 

In the equation of our electrical system, these elements 
will lead to terms analogous to the non-linear term in the 
Duffing Eq. ©. 

Fig. [4] shows the result of the complete transforma- 
tion of the initial circuit, in which the linear part of 
the circuit is now described by the admittance Y(u>) = 
Y e (u) + Zj 1 ^) in the Norton case, or by the impedance 
Z(iS) = Z e (uj) + Zj(iS) in the Thevenin case. 

Much qualitative insight can be gained from Fig. 
Indeed, it is clear that, in order to have the PNL par- 
ticipate and induce a significant nonlinear behavior, one 
needs the current going through it to be large. For this to 
occur, the linear admittance Y(uj) must be small. This 



will occur in the vicinity of a resonance frequency of the 
linear part of the circuit, namely when Im[V(|aj|)] = 0. 
The same reasoning applies for Thevenin representation 
(see Fig. HJj) where the voltage developed across the SNL 
must be large, i.e. in the vicinity of a resonance where 
Im [Z{u)\ = holds. 

B. The parallel and series LCR biasing circuits 

We now consider the simplest cases where the envi- 
ronment of the non-linear element - described by cither 
Y (id) or Z (ui) - is either the parallel or the series LCR 
circuit (see Fig. [S]). These circuits are defined by the 
inductances, capacitances and resistances L p , C p , R p in 
the parallel case and L s , C s , R s in the series case. From 
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TABLE I: Table of correspondance between the parameters 
of the various realizations of the generalized Duffing equation 
(see Eq. ((6])) considered in this article, i.e. the mechanical 
oscillator and the electrical oscillators based on a Josephson 
junction biased either through a parallel or a series LCR cir- 
cuit. The notion of participation ratio p has meaning only in 
the case of the Josephson junction circuits in which the cur- 
rent through the junction cannot exceed its critical current. 

the discussion in the last section, L p and L s include con- 
tribution of the Josephson effective inductance: for the 
parallel case, L~ 1 = L~l + Lj 1 and for the series case, 
L s = L se + Lj, where L pe and L se are the inductances 
of the embedding circuit of the physical junction. 

For the parallel case, the application of Kirchhoff 's law 
to currents in all the branches leads to the equation of 
motion 

C^o'S + ^S+^-sfl-^S 2 ) =/ rf cos^ (17) 
Up Lp \ bLj J 

which is a strict analog of Eq. ([2]). Here we have in- 
troduced the so-called gauge-invariant phase difference 
5 = $,/ (t) /ifiQ. In Table U we show the correspondance 
between the mechanical system and this parallel LCR 
system. The associated critical parameters are given in 
Table [Til Note that in this table the critical coordinates 
are referred to using the current I c = S c Iq. To a good ap- 
proximation, this is the amplitude of the current through 
the junction at the critical point. 

We now identify an important parameter which we call 
the parallel participation ratio 

p p = Lp/Lj (18) 

Together with the quality factor Q, it determines the 
ratio between the current I pc at the dynamical critical 
point and the maximum Josephson supercurrent Jo (see 
Table HI)) . The participation ratio measures the strength 
of the nonlinearity: a small participation ratio is as- 
sociated with a weak nonlinear term when d ~ 1 (see 
Eq. CLID). 

We now turn to the case of the series LCR circuit 
shown in Fig. Here, summing all the voltages across 
the elements of the circuit, we arrive at another equation 
of motion given by 

Lv^ + ljJI +Rq + Cq = V d cosvt (19) 



where q is the total charge having flown through the 
junction up to time t. The cubic nonlinearity of this equa- 
tion does not affect the restoring force like in the Duffing 
equation but the mass of the particle. Nevertheless, it 
is easy to show the same Duffing oscillator dynamics is 
recovered at low drives within the single harmonic ap- 
proximation (see Table |T|. Like for the parallel case, the 
critical amplitude in Table [TT] is given in a terms of a cur- 
rent I sc = u) c q c , which also here corresponds to the am- 
plitude of the current through the junction at the critical 
point. 

Again, we define the series participation ratio 

Ps = Lj/L (20) 

which plays the exact same role as p p for the series 
case. 

That the parallel and series case can be mapped into 
one another in the weak non-linear regime (see Table [TJ> 
is not an accident. Quite generally, one can show that 
any linear oscillator equation to which is added a cubic 
nonlinearity in any combination of x and its derivatives 
(i, x, etc.) will lead to, in the weak non-linear regime, 
the same dynamics as that of Eq. © . 

The striking conclusion of this section on a lumped 
element resonator is that even when the junction has a 
weak participation ratio, its non-linearity is not really 
"diluted". It will still display a bifurcation which can 
be employed for amplification, provided that the control 
of the amplitude of the oscillatory drive meets a corre- 
sponding increase in precision. 



C. General biasing circuit involving a resonator 

At the microwave frequencies where we wish to work, 
it is difficult to implement a pure LCR circuit without 
substantial parasitic elements. In practice, it will be eas- 
ier to implement a distributed element resonator built 
with section of transmission line o 19 ' 20 . However, as we 
are going to demonstrate, the conclusions of the last sec- 
tion are robust provided that the quality factor of the 
resonator is chosen adequately, which is easily achievable 
with on-chip superconducting thin film coplanar waveg- 
uides. We therefore now consider the general case of an 
arbitrary impedance Z (co) and use the Thevenin repre- 
sentation (see Fig. 0Jd). The extension of our results to 
the case of an arbitrary admittance in the Norton repre- 
sentation will be straightforward, using the set of simul- 
taneous transformation: Z(u) Y(w), Zj (lj) — ► Yj(u>), 
SNL^PNL. 

We start by writing the equations of motion using 
Kirchoff's laws and looking for a solution of the form 
I(t) = i/(£)e 1 "' + c.c By retaining only the drive har- 
monic terms and staying in the limit of small oscillation 
amplitude where we can take the nonlinearity into ac- 
count at the lowest order, one finds an equation of the 
form ©: 
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parameters at critical point 



response at critical point 



frequency 



drive amplitude 



oscillation amplitude 



mechanical 

parallel J+LCR 
series J+LCR 
parallel Y(oj) 
series Z(uj) 



s 2Q [ft - lj = -V3 
= 2Q - l) = -a/3 



£[M = -^(implicit) 



>[V(«c)] 

■[g("c)] 

Re[Z(u c )] 



Im[z( "-= )] = -a/3 (implicit) 



I d 

■*■ e 



3/2 

/ R.o[y(^)] \ 3/2 j 

5374 ^ |y>( Wc )| y j q 

8 / Rc[Z(^ e )] \ 3 / 2 
I \Zj(u, c )\ 



W4 
8 



f ILL 



/ c = <5 C J = pTTv/ ^¥ J o 



i w c g c = 



1 L r 



W c <£o 



j _ _J_ / R.c[y(^ c )] r 

Jc - 3 i/4y |yj(^ e )| i o 

J _ 4 / R.c[Z(^ c )] r 



TABLE II: Values of system parameters (frequency, drive amplitude) and response (amplitude of oscillation) at the critical 
point for the various mechanical and electrical realizations of the Duffing system considered in this article. In the two boxes 
marked "implicit", the critical frequency can be obtained only by the solving the given equation. 



n Z (n) / (n) (t) = I z 



-Zj 



J(t) 



I(t)-V d (21) 



^Here , Z = Z(uj = u d ), Zj = iuj d Lj, = 
d \<*>=Ud are t ne derivatives of Z(u>) taken at the drive 



d"/(t) 

dt" ' 



One can check that when 



frequency and j' n ) = 
I/Iq <C y/\Z/Zj\, i.e. when nonlinearity can be ne- 
glected, the equation describes the complete response of 
the system. In particular, the steady state oscillation 
amplitude is given by / = Vd/Z(cu). The zeros of Z(u) 
are thus the complex resonant frequencies corresponding 
to the normal modes of excitations of the linear system 
(A well-built resonator has only a sparse set of zeros reg- 
ularly distributed along the real frequency axis). The 
time derivative term on the left-hand side of Eq. (|2"Tj) 
accounts for the transient dynamics. 

The regime of interest for amplification is when 

I/Iq approaches y f \Z/Zj\ for frequencies in the vicin- 
ity of a zero of Z(ui). We then find an equation bearing 
a close resemblance with Eq. ([6]) for the first harmonic 
of the Duffing oscillator amplitude. In the following, we 
will neglect all higher derivatives Z ljl>1 \ 

It is important to note that while Eq. ^ yields the bi- 
furcation amplitude and drive as simple functions of the 
detuning frequency f2, equation (f2Tj) only leads to the am- 
plitude and drives Vd at the bifurcation points as a 
function of the complex quantity Z and Z' , which them- 
selves have complicated dependence on u> d - Although 
the expressions for generalized quantities such as /b(O), 
f§(Q) and fms(Q) are tractable, they provide little in- 
sight. We therefore refrain from presenting them in this 
paper. 

Nevertheless, we can still identify conditions that Z 
and Z' need to satisfy in order for the bifurcation real 
estate to exist. These conditions are obtained from Table 



Firstly, a critical drive frequency w c must exist 
such that the impedance Z(lo c ) satisfies the condition 
Re[Z(w c )] = — \/Slm[Z (u c )]. A way to meet this condi- 
tion is for the embedding impedance to have at least one 
resonant frequency ujq where Im[Z (wq)] = 0. Since for a 



passive circui 



-21 



dim [Z (lu) 



= Im [Z' (w )] > 



Re[Z(w )] > 
a critical frequency can be found at 



w 



Re[Z(w )] 



V3tm [Z 1 (w )] + Re \Z> (w )] 



assuming the impedance is a smooth function of fre- 
quency near the resonant frequency. In Fig. [5] we show 
how an exact solution for the frequency w c can be ob- 
tained graphically from a measurement of Im[Z (w)] and 
Re[Z(w)] in the neighborhood of ujq. The graphical con- 
struction shows that the critical drive frequency needs to 
be located on the low-frequency flank of a resonance of 
the bias circuit, as in the case of the simple LCR circuit. 

A second condition has to be fulfilled, however: the last 
line of Table HTl indicates that, at the critical frequency, 
the real part Re[Z(w c )] of the impedance must be much 
smaller than the effective impedance ui c (fo/Io = Ljlo c of 
the junction. The ratio of the RF critical current I c to 
the DC critical current Iq must be much smaller than 
one in order to fulfill the condition of weak non-linearity. 
The smaller the Rje[Z(ui c )]/\Zj(u) c )\ ratio, the larger be- 
comes the "real estate" available for bistability in the 
stability diagram. This last consideration sets the maxi- 
mum amount of dissipation which can be tolerated in the 
embedding impedance Z e (u>) . 

It is useful to translate this last condition in a more 
practical language. If we introduce the generalized qual- 
ity factor Q (w) of the environmental impedance: 
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Re[Z((u )] 






- Re[Z(ra)] 
■--j=Im[Z((»)] 

(0% 

I. 



FIG. 6: Sketch of the graphical determination of whether a 
circuit characterized by impedance Z(u) (see text) will al- 
low a Josephson junction with critical current To to bifur- 
cate. The procedure consists in i) plotting simultaneously 
Re[Z(uj)], — -^Im[Z(io)] and Zj(ui) = uj(po/Io, h) finding the 
critical frequency ui c (if it exists) at the eventual intersec- 
tion Re[zT(w)] = -~^=Im[Z(uj)], and iii) verifying that the 
condition for weak nonlinearity is satisfied at ui c , namely, 
Zj(u c ) > Re[Z(uj c )]. 



Re[Z env (u)] 



wIm [Z'env Ml 

wc find that the condition on I r can be rewritten as 



3V4V QW? 



< 1 



Thus, the higher the product of Q and p, the greater 
real estate for amplification. 



V. FURTHER DISCUSSION 

The phenomenon of bistability has been recently 
observed in an electrical circuit biasing a Joseph- 
son junctionals. Its operation as an amplifier has 
been studied^ and used to measure a superconducting 
qubrt^i 2 ^. The amplification near a bifurcation point is 
limited by a stochastic dynamical escape process similar 
to the one observed in DC biased JJ. This process is well 
studied theoretically in the classical regime where the 
fluctuations stem from thermal nois o 13 i 24 . The quantum 
regime is not as well understood but important results 
can be found i n 25 ' 26 . The circuit of Ref^ is based on 
the simplest realization of such an RF biased junction, 
namely, a parallel LCR circuit. The resistance was pro- 
vided by the 50 RF environment, the capacitance was 
constructed lithographically on chip, and the inductance 
came entirely from the junction itself. 



Another approach is to use continuous circuit elements 
instead of discreet ones such as L's, R's and C's. A res- 
onator structure with low internal loss can be used to bias 
the junction. Examples of these can range from a waveg- 
uide resonator, a lithographic resonator (see for e.g^-), 
or any other resonator geometry. Any harmonic of such 
structure can be used to bias the J J. The key lies in find- 
ing the correct way to couple the J J. Let us discuss briefly 
a particular implementation, namely, using a coplanar 
waveguide cavity resonator such as the ones recently used 
to observe cQED in a solid state system 2 ^. There, a junc- 
tion can be placed at the center of the transmission line, 
interrupting its central conductor. The circuit analysis 
shows that this corresponds to the case discussed above 
of a J J biased by a series LCR circuit. We have success- 
fully fabricated and tested such devices and will discuss 
the results in a later paper. 



VI. CONCLUDING SUMMARY 

In this paper we have presented and discussed the idea 
of harnessing, for amplification purposes, the bifurcation 
phenomenon of a driven Josephson junction embedded 
in a resonating circuit. We have first reviewed the ba- 
sic mathematics associated with the general saddle-node 
bifurcation phenomenon using the simplest example of 
the mechanical Duffing oscillator. We have then applied 
the corresponding formalism to a simple RF driven LCR 
electrical oscillator incorporating a Josephson junction 
and found the conditions for observing bistability and a 
diverging susceptibility. There are two relevant param- 
eters: the quality factor Q and the participation ratio. 
The latter quantity is necessarily smaller than unity and 
is a measure of the strength of the nonlinearity provided 
by the J J, relative to the embedding circuit. We have 
found that the product of these two parameters need to 
be much greater than one in order to provide a conve- 
nient operating conditions for amplification. We have 
then shown that the same condition applies to a more 
general resonating structure, with a simple adaptation of 
the definition of the quality factor and participation ratio 
to the vicinity of a resonance. Our prescription involves 
a relatively precise analysis of the microwave circuit of 
the resonator but can easily be followed with a varieties 
of geometries since microwave resonators can easily be 
made to have large quality factors. Later generations of 
amplifiers or detectors based on the bifurcation of super- 
conducting resonators should therefore benefit from the 
analysis presented in this article. 
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APPENDIX A: WEAK AND STRONG 
NON-LINEAR REGIMES FOR A NON-LINEAR 
OSCILLATOR 

In this appendix we would like to examine the behavior 
of a driven nonlinear oscillator beyond the simple frame- 
work developed in the main body of the paper which 
considers that the circuit only responds at the frequency 
of the driving force. We investigate the limits in which 
this hypothesis is violated, specializing to the realistic 
case of a Josephson junction. Our analysis relies heavily 
on the work developed i n 10 ! 11 ' 15 ! 16 ' 17 . We will consider 
the generalization of Eq. ([2]) to an arbitrary nonlinear 
potential. To simplify the notations and to avoid confu- 
sions with the main text, we use X — ► y, m — ► 1, replace 
LOd — > w and drop the noise term Fj^(t) from Eq. 
The purely non-linear force is written as N(y) and we 
thus start with 



y + 7y + uly + N(y) = f cos uit 



(Al) 



We will restrict our discussion to a nonlinear force 
N(y) which is an antisymmetric function of y as this 
is the case relevant to the nonlinearity provided by a 
Josephson junction. It will also simplify the discussion. 

A solution to Eq. (|A1|) can be found in the following 
form: 



y(t)=y° d %n) + 5y(t,n) 



(A2) 



where function y odd (t, n) is a sum of a finite number n 
of odd harmonics 



odd 



1 

k=l 



(A3) 



while 5y(t, n) is a correction to the expansion. Note 
that n is introduced here as an integer-valued parameter. 

By assuming Sy(t,n) -C y odd (t,n) we will conduct a 
linear analysis of this correction. Keeping Q = loq/^I 
constant (see Eq. (IA1[1 ). we find two different behaviors 
for Sy(t,n) depending on the parameters {oj,f} (their 
values are discussed later below) : 

1. Sy(t,n) is bounded in time and contains only odd 
harmonics of u) starting with 2n + 1. It follows that 
Eq. (|A3p is a good approximation for the stationary so- 
lution of Eq. (|Alj) . The choice of n is determined by 
the required accuracy of the solution. The amplitudes of 
the stationary oscillations y^ are functions of drive fre- 
quency lo and amplitude / given by a system of nonlinear 
algebraic equations in terms of N(y). 

2. 5y(t,n) is unbounded in time. This would lead to 
the breakdown of the validity of our linear analysis, or 
is at least an indication that the form for the solution 
given in Eq. (1A3|) has to be modified. The instability of 
5y(t, n) can be of two types: 



2a. 5y(t,n) ~ e At e «w(2&— 1)*^ wnere k < n is an integer 
and A > is a Lyapunov exponent. The instability of 
this type corresponds to the switching of the oscillation 
amplitude from one stationary state to another at fre- 
quency (2k — l)u>. This instability is of major interest to 
us as a resource for amplification purposes. In particu- 
lar, for k = 1, a switching between two stable states can 
occur. 

2b. Sy(t,n) ~ e At e i2fco>^ k < n is an integer. This 
is a different instability phenomenon because the solu- 
tion contains growing even harmonics which breaks the 
symmetry of the nonlinearity N(—y) = —N(y). It was 
shown to be a precursor of chaotic behavior of a nonlinear 
oscillator, at least for the Duffing case, where N(y) cx y 3 . 

Let us fix n and write down the system of nonlinear al- 
gebraic equations that defines yk, k = 1, 2, n by using 
the harmonic balance method: 



y 2fe _! [1 - (2k - lYu; 2 + i(2k - 1V 7 J + 

/ (A4) 

where <5i,2fc-i is the Kronecker delta function. The 
complex functions N2 m —i(zt, Z3, z-in-i) are defined as: 



A^ TO _l(zi,..,22n-l) = 

/ 2tt 



(A5) 



The solution of Eq. (j A4|) gives the amplitudes 
2/2fe-i(/j as a function of drive amplitude and fre- 
quency. Note that only the k = 1 harmonic is driven 
directly by an external force, while the higher harmonics 
feel the drive via the nonlinearity. 

The correction Sy(t, n) is obtained by subtracting the 
solutions to y odd (t, n) (see Eq. (|A4[1 ) from the solution 
of the initial equation (Eq. (IA1|) ) and keeping the linear 
terms in Sy. This leads to: 



Sy(t, n) + jSy(t, n) + u%6y(t, n)+ 

N'(y° dd (t,n))6y(t,n) = h(y° dd (t,n)) 



(A6) 



Here N'(y) = ™ 

One can see that the correction Sy(t, n) obeys 
the equation of a harmonic oscillator driven with 
a force h(y odd (t,n)) and parametrically driven with 
N(y odd (t,n)). The force h(y odd (t,n)) is defined by 



h(y° dd (t,n))= J2 N? k -i(yi,..,mn-i)e i{2k - 1)u,T + c.c. 

fc=n+l 

(A7) 

where c.c. is the complex conjugate. Note that 
h(y odd (t,n)) contains only odd harmonics starting from 



11 



2n + 1. If the oscillator is parametrically stable, the cor- 
rection Sy(t, n) will only contain oscillations at higher 
odd harmonics. They can be taken into account in prin- 
ciple by increasing the number of odd harmonics n in 
Eq. (|A3[) . This means that we can ignore the drive term 
for the analysis of parametric instabilities. The stability 
analysis will now be reduced to the analysis of the para- 
metrically driven harmonic oscillator which is very well 
understood. 

Because the function N'(y) is even N'(y odd (T,n)) will 
contain only even harmonics of u). That is 

oo 

N'(y° M (T,n)) = -^7V'L(yi,...,y 2 „-i)e^ t + C .c. 

(A8) 



fc=0 



with Fourier amplitudes N' 2k given by 

^'L( z b -,Z2n-l) 



z 2 fc-ie i(2fc - 1)fl +c.c. e- i2me 



de (A9) 



Now the equation for the parametric driven correction 
Sy can be rewritten in the form 



2 i TV' 1 



. 



0J* W 

where N' is evaluated at y\ . We get 
1 



Hi + -^Sy + (a + (3 cos 2t')Sy = 



(A12) 



(A13) 



this is known as Mathieu's equation with damping^. 
The main instability region corresponds to order one 
parametric resonance. It can be intuitively understood 
as a result of efficient pumping of an oscillator at the 
frequency y/a ~ 1. This instability leads to growing os- 
cillations of Sy ~ e A * e l * = e Xr e tuJt . Importantly, there 
are at the same frequency as our odd anzatz y odd (t,l). 
This can be incorporated into the anzatz by replacing 
yi — ► yi(t). This was done in the main body of the pa- 
per and was shown to lead to hysteresis and bistability. 
In terms of a(u>,y~i) and f3(ui,yi) the unstable (bistable) 
region is given by 



a{u,yi) = 1 ± 



Q 



(A14) 



32 



Sy(t, n) + jSy(t, n) + Sy(t, n) x 

r oo \ 

1 + #'o + £ (^2*(Wi. ->ihn-i)e aiut + cc.) = 

v k=l J 

(A10) 



Eq. (|A10[) is known in the literature on differential 
equations as Hill's equation with linear damping. Meth- 
ods to investigate its stability diagram both analytically 
and numerically can be found inii for example. 

The goal of our discussion can be achieved by consid- 
ering the simplest form of (|A10|) . We will take n = 1 
corresponding to the single mode solution y odd (t, 1) = 
yi(t) = \{y\e lut + cc.) and truncate the series (|A8[) to 
only the first and second terms. This corresponds to a 
DC shift in the linear oscillation frequency and to the 2u> 

parametric drive with amplitude N' 2 {yi). This approx- 
imation is rich enough to understand the bistability of 
the first harmonic response of a driven nonlinear oscilla- 
tor as discussed in the main body of the paper. It also 
shows the roads that lead to the breakdown of the simple 
picture of bistability. These simplifications lead to 



Sy + jSy + 5y[ul 



N'l(yi) 



\N'l(y 1 )\co S (2ut + Arg[N>l})} = 



(All) 



where we have simplified our notation by using 
8y(t, 1) = Sy. To reach the canonical form let us shift 

and rescale time by defining 2t' = 2tot + Arg[N' 2 ] and 
introduce two parameters of central importance: 



The two nearest instabilities of the "wrong" type cor- 
respond to a second order parametric resonance when 
i/a ps 2 and to the phenomenon of the drive-mediated 
negative restoring force when a < 0. These two instabili- 
ties contain growing double- frequency and zero- frequency 
components in Sy(t), respectively and they break the 
symmetry of the non-driven problem. The locus of these 
transitions, with an accuracy of order o(\/3(u>, yi)| 2 ) is 
given respectively by 



and 



1 



r, ^2 



<*(v,yi) = --/3(w,j/i) 



(A16) 



For the Duffing potential, where N(y) — w 2 ?/ 3 we 
found that a and /3 have particularly simple form: a — 

We now arrive at the main result of this Appendix: 
provided that \f3 D \ = ^||y 2 < 1 and 1/2 < w = uj/aj < 
1 the symmetry breaking unstable regimes are inaccessi- 
ble since neither of Eqs. (|A15|) . (|A16|) has a solution. The 
analysis conducted in this Appendix provides the ground 
for discussing the nonlinear oscillator only in terms of a 
single frequency oscillating response. We saw in the main 
body of the paper how this leads to a simple picture of 
bistability. The exact boundary at which this picture 
breaks down requires a sophisticated analysis. However 
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the simple analysis reveals that the instability regions 
are well separated from our region of interest for the 
relevant experimental parameters (1 — w ~ 1/Q <C 1, 
\yi\ ~ ^ !)■ More importantly, it shows that 

the breakdown of the simple bifurcation picture is Q- 
independent in the limit of high Q and the size of the 
accessible bistability region can therefore be under con- 
trol. 

In conclusion, basing ourselves on refer- 
ence a 10 ' 11 ^ 15 ' 16 ! 17 we have explained that instabilities of 
the steady state response of a Josephson junction non- 



linear oscillator can be viewed as parametric resonances. 
This mapping allows to classify the different instabilities 
and separate the ones of interest for bifurcation amplifi- 
cation from the unwanted ones. A "sound" bifurcation 
is garanteed when i) the quality factor is much greater 
than unity, ii) the relative detuning is of the order of the 
inverse of the quality factor and iii) the dimensionless 
oscillation amplitude is of order of the inverse of the 
square of the quality factor. 
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